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Abstract 

This paper proposes an empirical test of financial contagion in European equity 
markets during the tumultuous period of 2008-2011. Our analysis shows that tradi- 
tional GARCH and Gaussian stochastic-volatility models arc unable to explain two key 
stylized features of global markets during presumptive contagion periods: shocks to 
aggregate market volatility can be sudden and explosive, and they are associated with 
specific directional biases in the cross-section of country-level returns. Our model re- 
pairs this deficit by assuming that the random shocks to volatility are heavy-tailed and 
correlated cross-sectionally, both with each other and with returns. The fundamental 
conclusion of our analysis is that great care is needed in modeling volatility if one wishes 
to characterize the relationship between volatility and contagion that is predicted by 
economic theory. 

In analyzing daily data, we find evidence for significant contagion effects during the 
major EU crisis periods of May 2010 and August 2011, where contagion is defined as 
excess correlation in the residuals from a factor model incorporating global and regional 
market risk factors. Some of this excess correlation can be explained by quantifying the 
impact of shocks to aggregate volatility in the cross-section of expected returns — but 
only, it turns out, if one is extremely careful in accounting for the explosive nature of 
these shocks. Wc show that global markets have time-varying cross-sectional sensitiv- 
ities to these shocks, and that high sensitivities strongly predict periods of financial 
crisis. Moreover, the pattern of temporal changes in correlation structure between 
volatility and returns is readily interpretable in terms of the major events of the periods 
in question. 
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1 Contagion and financial crisis 

Contagion refers to the idea that asset returns in inter-related financial markets sometimes 
exhibit anomalous patterns of correlation. Yet these anomalies are difficult even to identify, 
much less characterize empirically, and their precise causes — whether financial, macroeco- 
nomic, or behavioral — arc hotly debated. 

In this paper, we study the period of presumptive financial contagion associated with 
the European sovereign-debt crisis of 2008-2011, which arose out of market fears that 
Greece would default on its sovereign debt. The crisis soon grew to include fears about 
the future solvency of other European nations with large debt loads. These fears, and the 
associated turmoil in global financial markets, led to the creation of the European Financial 
Stability Fund (EFSF) in May of 2010, which was intended to faciliate low-cost loans for 
other struggling EU members, including Portugal and Ireland. Yet uncertainty continued to 
cloud the global macroeconomic outlook throughout 2011 and early 2012, a period during 
which the likely impacts of Eurozone contagion on financial markets, monetary and fiscal 
policy, and trade links were far from known. 

During this period of instability, stock returns in both European and US markets exhib- 
ited explosive daily movements that are very hard to explain using traditional asset-pricing 
models. To understand the magnitudes of the effects that have to be modeled, consider 
that the German DAX index fell from a high of 7402 at the beginning of July 2011 to 5216 
by 4 October, resulting in a -30% return. Over the same period the UK FTSEIOO had a 
fall of similar magnitude. Prom a volatility perspective, the S&P 500 volatility index (or 
VIX) dramatically increased from 15% on July 1, 2011 to 45.5% by 4 October, staying 
persistently high for at least a week thereafter. 

The sheer length of the European debt crisis provides a natural testing ground for 
models of contagion. We take an asset-pricing perspective on this problem, and propose a 
multifactor, explosive-stochastic-volatility modeling framework. This approach generalizes 
existing models by allowing for the possibility of explosive, mutually exciting shocks to 
volatility that play a role in the cross-section of expected returns. Our work therefore builds 
on the long line of research into stochastic volatility and inference for Bayesian state-space 



2 



(1993 


) Kim et al. 


(1998 


), 


Aguilar 



and West (2000), Liu and West (2001), Jacquier et al. (2002), and Chib et al. (2005). Other 



important references on Bayesian computation for non-linear state-space models include 



Carvalho et al. (2010a), Niemi and West (2010), and Lopes and Tsay (2011), among others. 



The closest related model is that of Lopes and Migon (2002), who also incorporate 



stochastic volatility into an asset-pricing model for contagion. The crucial difference in 
our model is that it can express all three ways in which market shocks tend to cluster 
during times of crisis: time-series clustering, where large shocks today correlate with large 
shocks tomorrow; cross-sectional clustering, where large shocks in one region correlate with 
large shocks in other regions; and directional clustering, where shocks to aggregate market 
volatility correlate with specific directional biases in contemporaneous country-level returns. 
These first two forms are, respectively, the "heat wave" and "meteor shower" metaphors 



of volatility (or jump) clustering proffered by Engle et al. (1990). (See also Eraker et al. 



2003: Ai't-Sahalia et ahj 2010) These forms of clustering are certainly not unique to con- 



tagion periods. Our approach differs from this line of thought, in that we emphasize the 
role of directional clustering, driven by aggregate volatility shocks, in explaining contagion 
episodes. To cite one example: during the European sovereign debt crisis of 2010, our anal- 
ysis finds that aggregate volatility shocks were systematically associated with large positive 
residual returns on the Spanish stock market, and large negative residual returns on the 
German market. Existing volatility models are agnostic with respect to the sign of these 
residuals, and provide no mechanism for aggregate volatility shocks to be associated with 
specific directional biases in country-level moves. Yet such biases show up clearly in the 
data. 

We analyze daily data on returns from country-level equity indices during the European 
sovereign debt crisis, and reach three main empirical findings. First, we present evidence 
that European equity markets exhibited significant excess correlation during the debt crisis 
of 2010, relative to an asset-pricing model that assumes regional market integration. Figure 
[T] gives a brief visual summary of this evidence. It depicts the trailing two-month residual 
correlation coefficients among Germany, Italy, and Spain between May 2008 and October 
2010. The residuals are from a two-factor model incorporating global and European market- 



Contagion and Volatility, May 2008 - October 2010 




Figure 1: Top: trailing two- month correlation coefficients for the residuals of Germany- 
Spain, Germany-Italy, and Italy-Spain from a two-factor model. Middle: our proposed risk 
factor that measures regional excess volatility. Bottom: estimate of global market volatility 
that arises from our explosive stochastic-volatility model. 
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risk factors, described in Section [5j In April and May of 2010 (shaded grey in the figure), 
there is a clear divergence from the historical norm, precisely coinciding with the Greek 
sovereign-debt crisis and associated bailout. 

Second, we find that part of this excess correlation can be attributed to the impact 
of mutually exciting volatility shocks in the cross-section of expected returns. We propose 
global and regional volatility-risk factors to quantify this impact, and show how to construct 
these volatility risk factors directly from the time series of market returns. This is similar 



to the work of Ang et al. (2006), but differs in two key respects: we consider returns at the 
level of national portfolios, rather than individual domestic equities; and we construct the 
volatility risk factors directly from the time series of market returns, rather than indirectly 
from option prices. 

Our results show that markets have time- varying sensitivities to these risk factors, which 
implies time- varying correlation between shocks to volatility and contemporaneous shocks to 
returns. Crucially, periods of increased sensitivity to these risk factors coincide with periods 
of excess correlation. It is only a slight oversimplification of our findings to say that volatility 
shocks provide little predictive value during quiet periods, but significant predictive value 
during times of market crisis, which are precisely the times when market shocks exhibit 
strong directional clustering. Figure [T] shows this correspondence: the regional volatility- 
risk factor accurately predicts periods during which correlations among the residuals for 
Germany, Italy, and Spain diverge significantly from their long-run average. The pattern of 
volatility factor loadings that we find during crisis periods, moreover, is readily interpretable. 
For example, during May 2010, the regional volatility factor partially accounts for two 
notable phenomena: Spanish, Italian, and Belgian markets rallying in anticipation of a 
Greek bailout; and German and British markets suffering in anticipation of footing the bill. 

Our third main finding concerns the contemporaneous relationship between volatility 
and expected returns. This finding is tangential to our main point, but is of independent 
interest in the finance literature. We are primarily interested in the role of volatility shocks 
in understanding contagion, rather than in estimating volatility per se. Nonetheless, in 
performing a formal statistical assessment of our volatility model, we find a significant 
negative relationship between daily volatility and expected returns on the U.S. market 



portfolio. Crucially, no such finding on tfie daily time scale can be supported using other, 
more traditional volatility models. This finding links our methodology with the longstanding 



debate on the relationship between risk and return (e.g. Glosten et al. , 1993). 

Our empirical approach relies upon a new method for measuring volatility shocks, which 
is a necessary step in constructing the volatility risk factors that we propose. As we will 
show, traditional volatility estimators cannot explain the large shocks to volatility seen 
during recent financial crises, and sometimes lag several days behind during times of acute 
turmoil. Many such large volatility spikes are evident in Figure [TJ This leaves investors 
in the unsatisfying position of accepting that asset prices undergo moves of 6, 8, or even 
10 standard deviations far more often than their models would predict. To correct this 
problem, we propose an explosive stochastic volatility (ESV) model, where the random 
shocks to aggregate volatility arise from a distribution whose tails are much heavier than 
Gaussian, and which can be fit straightforwardly using particle learning. This model does 
not require implausibly large shocks in order to generate observed movements in asset 
returns. 

The rest of the paper proceeds as follows. In Section 2, we describe the role of volatility 
in understanding presumptive periods of financial contagion. In this context we review 
several, related proposals that motivate our decision to focus on how volatility shocks enter 
the cross-section of expected returns, providing important context for our proposed asset- 
pricing models. In Section 3, we explain the most general version of our mutually exciting 
stochastic volatility model, wherein the correlation between volatility shocks and return 
shocks are explicitly quantified. We also describe a handful of empirical simplifications 
that we make in order to fit the model to country-level European returns. Our estimation 
strategy then unfolds in two stages. In Section 4, we estimate both the expected returns and 
the daily volatility states for both the U.S. and European markets using particle learning. 
We also describe the results of an extensive model assessment, wherein we benchmark our 
volatility model against an alternative model based on the GARCH framework. In Section 
5, we use the filtered mean and volatility states from Section 4 to fit the encompassing 
model described in Section 3. Sectino 6 concludes with some final remarks about how our 
findings fit into the existing literature. 
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2 The role of volatility in understanding contagion 

Existing economic theory provides several reasons why the study of market volatiUty — the 
primary modehng focus of our paper — would play a major role in understanding contagion. 
We divide these into two broad categories: economic (or fundamentals-based) reasons; and 
statistical reasons. 

Many previous authors have provided economic interpretations of why volatility shocks 
would predict variation in asset returns during crisis periods, even when other risk factors 
have been priced. One explanation is that volatility shocks are simply a proxy for factors 



that can be loosely grouped into the category of investor behavior (e.g. Dornbusch et al. 



2000 ) . A second possible interpretation can be found in the literature on the portfolio flows 



of large international investors. Froot et al. (2001) study the relationship between local 



market returns and capital inflows from large institutional investors domiciled outside the 
local market. They find evidence that these inflows are typically associated with larger 
expected returns, but that transitory inflows predict lower returns. To the extent that large 
regional volatility shocks predict increasingly transitory capital flows, they may also predict 
returns. 

A third possible interpretation is that volatility shocks are a proxy for changes in the 
informational efficiency of equity markets. The resulting effects may take several forms. 



Calvo and Mendoza (2000) observe that, due to the potentially high cost of procuring and 
assimilating country-specific information, investors may rationally engage in herding be- 
havior due to mutual dependence upon a relatively small cohort of information conduits. 
Additionally, standard models of market microstructure suggest that investors with infor- 
mational advantages will accrue their desired market positions slowly, in order to minimize 
transaction costs. Therefore, the portfolio flows of informed investors will be more highly 
autocorrelated than those flows corresponding to uninformed transactions, and therefore 
more strongly predictive of future returns. If informational asymmetries are high, flows 
involving investors with an informational advantage will constitute a larger fraction of total 
flows. Under this assumption, if volatility shocks correspond both to information asymme- 
tries and to increased trading volume, they will also enter the cross-section of returns. 
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Moreover, stochastic volatility offers perhaps the best modeling framework for under- 
standing contagion when accurate data on macroeconomic fundamentals is lacking. The 
absence of data on fundamentals is especially acute when studying returns on a daily scale, 
rather than a scale of months or weeks. Working with daily data poses unique challenges: 
the potentially large size of daily volatility moves, the poor signal-to-noise ratio, and the 
fact that different markets operate on different trading calendars. The models we use in 
this paper directly address the first two issues. To address the third, we use daily data 
on U.S. -based exchange-traded funds (ETFs) corresponding to major European equity in- 
dices, rather than the indices themselves. The use of market proxies constructed from ETFs 
(which are highly liquid and free of any apparent arbitrage on a daily time scale) eliminates 
the complications due to differences in the trading calendar across different nations. 

Using data on coarser time scales potentially allows important information such as trade 



variables and interest rates to be incorporated into the asset-pricing models, as in Bekaert 



et al. (2005). This allows for a more precise economic interpretation of any contagion effects 



that might be found. Yet low-frequency data imposes an even larger disadvantage, given 
the poor statistical precision it allows in separating contemporaneous from time-lagged 
covariation. Since trade conditions and other macroeconomic variables do not change fast 
enough to drive much of the variation in daily returns, we are left with volatility shocks as 
the most salient estimable quantity to use in our pricing model. 

Finally, there are purely statistical issues associated with identifying contagion. The 
most difficult of these are related to the issue of time-varying volatility. Suppose that one 
has specified an asset-pricing model that captures all relevant shared risk factors for a group 
of assets, such that the idiosyncratic components of variation are uncorrelated for all assets. 
Even under this ideal scenario, a deceptive illusion of excess correlation may arise if the 
presence of time- varying volatility is ignored. 

Such an effect can be seen in the context of an international version of the capital asset 
pricing model (CAPM) with time-varying volatility. Suppose that yn is the excess return 
of the national equity index for nation i at time t; that fj,t is the conditional expected excess 
return of the world market, given information available at time t — 1; that xt is the realized 
excess return on the world market; and that xt = fj,t + et, where var(et) = ct|. If markets 



are globally integrated in accord with the world CAPM, then 



Hy^t) = Pi^t , (1) 

where the residuals eu = {yu — l^ifJ-t) have variance af^. 

Even in this simple setting, ignoring the effect of time-varying volatility can lead to 
spurious findings of excess correlation. This happens because cross-market correlations are 
conditional upon the overall market volatility , since the market shock et affects the return 
for every asset. Hence when market volatility is high, estimates of these correlations are 



biased upward (see, e.g. Forbes and Rigobon 2002). A number of authors have attempted 



to correct for this in studying previous crises previous crises. For example, [Bekaert et al. 



( 2005 ) explore the importance of segmentation versus integration in explaining contagion, 
and consider the possibility that increased volatility of regional risk factors can account for 
excess correlation during the 1997 Asian crisis and the 1994 Mexican crisis. 

The problem is further exacerbated by an errors- in- variables effect. Suppose that one 
were to fit the CAPM using the realized market returns xt . Then the least-squares residuals 
corresponding to assets i and j will be correlated, even if all variance terms are known, and 
even if the real e^s are cross-sectionally independent. This excess correlation, moreover, 
will be higher for those residuals corresponding to times of market stress, when is large. 
To see this informally, note that 



cor(yit, yjt) = cor{/3j(xj - et) + eu, Pj{xt - et) + ejt} 
PA 

which is large when erf is large relative to the idiosyncratic variances af^. This excess residual 
correlation results from two well-known facts: that the least-squares estimators for /3j and 

are biased downwards when the expected market return is observed with error; and that 
this bias is more severe when cj^ is large. 

Therefore, even asset-price movements that arise from the most basic form of the CAPM 
can look suspiciously like excess correlation, unless both the errors-in-variables problem 
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and the time- varying volatility problem are handled appropriately. This kind of correlation 
is best thought of as correlation arising from integrated markets, and thus is not really 
contagion according to the asset-pricing definition used in this paper. 



3 Mutually exciting stochastic volatility 

3.1 Motivation 

We therefore adopt the goal of quantifying the effect of volatility on market returns during 
presumptive contagion periods, remaining agnostic on the causal interpretation of these 
effects. We focus most specifically on the phenomenon of "directional clustering" referred 
to in the introduction: namely, that a large shock to aggregate volatility predicts specific 
biases in the signs of country-level returns. Figure [T] hints at precisely this effect, when, for 
example, the May 2010 spike in Eurozone volatility was associated with large, historically 
aberrant negative correlations between Germany and Spain. 

Existing models do not explicitly address this effect, whereby epochs of high volatility 
are associated with specific directional biases in country-level asset returns. For example, a 
traditional multivariate stochastic- volatility model asserts that a large positive residual for 
the Italian stock market during May of 2010 was equally likely to be associated with a large 
positive residual for Germany as it was a large negative residual. Figure [T] puts the lie to 
this assumption. In contrast, our goal is to characterize these directional biases explicitly, 
and to understand the role they play in presumptive contagion episodes. 

3.2 Modeling framework 

Let Uit represent the excess returns for asset i during period t. A general, dynamic multi- 
factor representation of expected returns would take the following form: 



Here fjt represents a shared risk factor, with fijt denoting the corresponding conditional 
mean of that factor. For example, the CAPM sets K = 1, ft as the market return, and fj,t 



K 




(2) 
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as the expected market return. 

Our model specification accounts for the conditional volatility of the risk factors them- 



selves, taking the square-root volatility model of Heston ( 1993 1 as a starting point. Through- 



out the following exposition, we assume that all volatility states are restricted so that the 
square-root volatility process never becomes negative. 

We describe the most general formulation of our model, before considering specific sim- 
plifications that are necessary in order to fit the model to data. In matrix form, we assume 
that 



ft = f^t + Ut^t 
yt = cxt + Bt{ft - fJ't) + Vtet 



(3) 
(4) 



where it = {fit, ■ ■ ■ , fxt)' is a vector of shared risk factors at time t; Bt is a matrix of 
factor loadings; = {yit, ■ ■ ■ iVpt)' is a p-vector of excess returns for individual national 
equity markets; and where Ut = diag(uit, . . . , UKt) and Vt = diag(vit, . . . , Vpt) are diagonal 
matrices of volatility states (i.e. conditional standard deviations). 

We assume a multivariate stochastic volatility model for the complete vector of volatility 
states (uf, Vt) = {uu, ... , UKt,vit, • ■ • , Vpt)': 





( 



\ 



UU 

^> 







' Vt 



+ 



(5) 



Vt 



where rj^ = (r/f, r]^)' is a vector of longitudinally independent innovations to volatility. See, 
for example, [Eraker et al] ( |2003| ) and |Chib et al.| ( |2005 ) . 

Finally, define i/?^ = (^(, et, rj^)' to be the stacked cross-sectional vector of innovations at 
time t in the factors, returns, and volatility states, respectively. Our model is completed by 
specifying a covariance structure for these innovations, which we assume to be multivariate 
Gaussian: 

V't~N(0,Et). 
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Suppressing the time index for ease of notation, partition S as 
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Some further identifying restrictions on are necessary, due to the presence of the 
diagonal scale factors Ut and Vt in the state evolution equations, and because the factor 
innovations also explicitly contribute to the observation equation Q, since it — fJ-t = Ut$,t. 
To address these, we assume that: 

1. E^^ and S^e are correlation matrices; and 

2. the factor innovations are cross-sectionally independent of the return innovations 
et and the volatility innovations 'q^, so that = Ti^u = '^s,v = 0. 

We are therefore left with a simplified model where 
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with inverse partitioned conformally as 
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Even with these identifying restrictions, the resulting model is still very flexible. Impor- 
tantly, it can accommodate all three forms of clustering described above. 

Cross-sectional: A large shock in one market can propagate to other markets in two ways. 
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First, volatility shocks are themselves cross-sectionally correlated at level of factors 
(Sum), returns {Ti^v), and factor-return interactions (Sui,). Second, time-lagged cross- 
sectional correlations can be amplified by the regression matrix D in the volatility- 
state evolution equation. 

Time-series: Equation ([s]) descibes a state-evolution equation wherein volatility states are 
auto- and cross-correlated over time. 

Directional: Unless both Ti^u and T,^^ are identically zero, shocks to volatility will be 
associated with specific biases in the signs of the return residuals (ej). 

The theory of the previous section predicts that periods of contagion would be associated 
with changes in and — that is, the correlation between shocks to volatility and shocks 
to returns. We are now in a position to formally investigate this prediction. 

3.3 Estimation strategy 

We view the general dynamic model we have just outlined as a useful conceptual framework 
for thinking about contagion. Nonetheless, it poses a number of difficult statistical issues 



associated with high-dimensional state-space models (see, e.g. Liu and West , 2001 Carvalho 



et al. 2010a). These include the unknown parameters in the state evolution equations; the 
large dimension of the full covariance matrix S^; the potentially heavy-tailed character of 
the volatility shocks; the errors-in-variables problem that arises because the conditional 
factor means are unobserved; and the need to specify a model for how parameters such 
as Bt change over time. 

Therefore, rather than attempting to solve all of these issues in the context of a single 
dynamic linear model, we make a number of specific simplifications for the purpose of fitting 
this model to data from the European sovereign debt crisis. These simplifications, and our 
two-stage estimation strategy, are guided by the desire to test the main prediction of the 
full model: namely, that innovations in the volatility sequence enter the cross-section of 
returns (yt), above and beyond their contribution to the volatility of the factors {ft)- 

To see this, suppose we are given the values of the volatility innovations r/j. Appealing to 
standard multivariate-normal theory, the conditional distribution of the return innovations 
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et IS 



where 



-1 



rrii 



\ 
( 



yl y 



(6) 



(7) 



Therefore we may write the conditional model for returns as 



yt = cct+ Btift - lit) + rr^r + + vt^t 



(8) 



where 



are matrices of regression coefficients, and where var(et) = H^. 

Observe that, in the conditional model of Equation (|8]), the volatility innovations behave 
like shared risk factors that enter the cross section of returns. We can therefore describe 
the model in either of two equivalent ways: 1) time- varying correlation between residuals 
in the returns and volatility evolution equations, described by blocks ^eu,t and ^ev,t of the 
overall precision matrix fit = S^^; or 2) time- varying sensitivities to volatility-risk factors, 
described by loadings Tf and FJ'. 

We work directly with this second (conditional) model, subject to several further simpli- 
fications. First, we assume a fixed-parameter, two-factor model that describes daily excess 
returns in terms of exposure to global-market and regional-market risk factors, along the 



lines of Bekaert et al. ( 2005 ) . This leads to the following marginal model for the returns 
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series: 



„ -^j_aUS(fUS ,US\ j_ oEU( fEU 



where yu is the excess return for index i on day t; ff^ (fJ-Y^) is the observed return (expected 
return) on the value-weighted US market portfoho; and ff^ (/^f ^) is the observed return 
(expected return) on the pan-European index from Morgan Stanley Capital International, 
which corresponds to a value- weighted portfolio drawn from the 16 largest equity markets 
in Europe. We use the US market return as a proxy for the global market. We also use the 
return Vanguard's European index exchange-traded fund (ticker VGK) as a proxy for the 
MSCI European index. Since this ETF is traded in the US market, it allows us to sidestep 
differences in trading calendars among the individual European markets themselves. 

Second, we assume that the full covariance matrix S (which describes all correlation 
among the shocks to factors, returns, and volatility states) is restricted in such a way that 
Qev = for all t. This imposes the requirement that shocks to country- level returns are 
conditionally independent of shocks to country-level volatility, given the shocks to factor 
volatility. We have therefore specified a particular form of Gaussian graphical model for S 



(see, e.g. Carvalho and West, 2007). Accordingly, rj"^ drops out of the conditional model Q 



for returns, with only rjf remaining (in addition to the factors themselves). 

Third, we assume that the innovations to country-level volatility and factor-level volatil- 
ity are conditionally independent, given all other contemporaneous parameters and states of 
the DLM. This encodes the assumption that the factor-level volatility innovations are suffi- 
cient to describe the correlation between volatility and returns, and that the country-level 
volatility shocks offer little extra information. 

Fourth, we assume that the factor volatilities (n|,?i2) for the US and EU markets, 
respectively, obey the following triangular system: 



Ut 



V 



+ 



V 











V 



H-l 



H-l 



+ 
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hi 
hi h2 



J 



5| 



(9) 



where 5j and 5^ are mutually independent. We interpret 5i as a shock to global market 
volatility, and 6^ as an independent shock to region-level volatility. The actual volatility in 



15 



the regional market-risk factor is then a Unear combination of 6} and 6^, since from ([& 



^ hi \ 
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This imphes that the lower-triangular matrix L is the Cholesky factor for the marginal 
covariance matrix of the factor innovations rj^. 

Finally, we assume that shocks to returns are conditionally independent of each other, 
given the factors, volatilities, and volatility innovations. Statistically speaking, this means 
that il.ee is assumed to be a diagonal matrix. 

We are therefore left with a simplified form of the full model wherein 
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This, in turn, implies that the conditional model for returns ([8]) reduces to 

= a + B{it - fit) + r^Vt + Vet , (10) 



where the et's are independent of 77" with variance 



var(ei) = 9. J = See - Se„S , 



These simplifications suggest that, if one knew the volatility states, it would be possible 
to estimate the time-varying pattern of correlation between volatility shocks and return 
shocks, strictly by estimating F". This leads us to adopt the following two-stage approach 
for estimating the model. 

1. Fit independent heavy-tailed stochastic volatility models to the US and EU markets, 
respectively, for the purpose of getting filtered estimates for the conditional means 
and volatilities of the factors. This is described in the next section. We use these 
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filtered states to regress the EU volatility innovations on the US volatility innovations 
in order to extract the two independent sets of volatility risk factors. 



2. Then estimate the conditional model in (10), fitting a traditional factor model using 



the volatility innovations as shared risk factors. By using 6j and 6^ rather than rjj 
and r/l, we eliminate collinearity among the predictors. 

4 Estimating time- varying market volatility 
4.1 The relationship between risk and return 

Our first, intermediate goal is to model the conditional mean and volatility of the daily 
return on the US and EU markets. This allows us to isolate the global and regional market 
shocks, and to extract a volatility risk factor for incorporating into the pricing model for 



the market returns of individual national equity indices (10). It is also essential if we are 
to correct for the errors-in-variables effect that arises because the factor means fi^ in this 
model are unknown. 

Extracting a volatility-risk factor requires specifying the the fundamental relationship 
between risk and expected return on the market. The market models we consider are 
simple and fairly standard, and are identical to those used in the literature on the equity 
risk premium. Yet as we will see, they are still tricky to estimate. 



Following Merton (1980), Glosten et al. (1993), and others, we consider three baseline 



models for the excess market return. The first two of these models also have special "no- 
intercept" cases, which we denote Models 1(b) and 2(b), respectively, giving five overall 
models. To avoid confusion with the notation already established, we phrase the issue 
somewhat generically, letting xt denote the market return and at the conditional standard 
deviation of the market return. The models we consider are: 

Model 1(a): xj = oq + aiaf + et , var(e() = 
Model 2(a): xt = ao + aiat + et , var(et) = at 
Model 3: xt = ao + et , var(et) = af 
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Model 3 hypothesizes that expected market returns are constant over thne, and are unrelated 
to the conditional volatility of the market. We also entertain Model 1(b), the special case 
of 1(a) where ao = 0. This approximates an equilibrium model in which the relative 
risk-aversion function for the representative investor is constant over long periods of time. 
Likewise, we entertain Model 2(b), the special case of 2(a) where ao = 0. This restriction 
implies that the market price of risk, E(xj)/c7(, is constant over long periods of time. 

Rational, risk-averse investors require higher expected returns to compensate them for 
taking on greater levels of risk. Yet this need not imply a positive relationship between risk 
and return over long time periods. For example, if investors systematically tend to save 
more of their wealth during risky times, then asset prices may be bid up considerably, and 
expected excess returns may go down as a result. For this and other reasons, the sign of 
ai may be either positive or negative without contradicting the assumption that investors 



are rational. See Glosten et al. (1993) for an overview of the literature on the relationship 
between market volatility and the equity risk premium. 

These models have a relatively long history in finance, where the need to account for 
time-varying volatility in characterizing the expected market return has long been recog- 



nized. For example, Merton (1980) considers Models 1(b), 2(b), and 3 for monthly data, 
using very simple estimates of monthly volatility. Merton concludes his study by observ- 
ing that "[the] most important direction is to develop accurate variance estimation models 



which take account of the errors in variance estimates." Glosten et al. (1993) do precisely 
this. Focusing on Model 1(a), they use an asymmetric GARCH model incorporating the 
risk-free interest rate and seasonal dummies to estimate monthly volatility. They find evi- 
dence that ai is significantly negative, implying a negative conditional correlation between 
returns and volatility. 

Many further elaborations of this basic framework are possible, but typically necessitate 



working on a coarser time scale. Bekaert et al. (2005), for example, incorporate many 



economic fundamentals into their model of the the expected market return. We do not 
incorporate these fundamentals into our models because they change very slowly on a daily 
time scale and will therefore drive very little of the variation. 

Each of the models above requires further assumptions about the behavior of the volatil- 
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ity term at- We assume a square-root stochastic volatility model with heavy-tailed shocks: 



(Tt = (Jo + (l)at-i + Xtzt (11) 
At ~ tt (12) 
Zt ~ N(0,r2), (13) 

with at truncated below at zero, and where is a positive-t random scale with degrees 
of freedom. 

This differs from a traditional Gaussian SV model, in that the innovation is itself a 
product of two terms: a normal shock (which is present in a normal SV model), and a 
half-t-distributed "local" scale factor (which isn't). This results in a horseshoe distribution 



(when V = 1) for the marginal of the shock XtZt (Carvalho et al. , 2010b), with tails far 
heavier than Gaussian. Formally, it corresponds to a discretized version of subordinated 
Brownian motion for the square-root volatility process, a construction first explored by 
S. Bochner in the 1950's. The increments of the random clock, or subordinator, are half-t 
distributed. 

The intuition behind our approach is that the daily volatility of asset returns is highly 
persistent during normal times, but on rare occasions will change drastically and rapidly. 
When these drastic moves occur, rational investors must quickly discount the past and 
revise their expectations about volatility. 

Fat-tailed versions of traditional stochastic volatility models have been proposed to un- 



derstand volatility filtering Jacquier et al. (2002) and currency returns in a multi-factor 



setting Aguilar and West (2000). Our approach here is different in that we explore the 
expected-return relationship with time- varying volatility. This allows us to provide a volatil- 
ity factor that also respects the fundamental risk-return trade-off inherent in an asset-pricing 
market model. We also fit the model with a fundamentally different purpose: as the first 
stage in our two-stage estimation procedure for quantifying the cross-sectional relationship 
between volatility shocks {Xt^t) and return shocks. Our empirical findings show strong ev- 
idence of fat-tails, and we are able to use this conclusion to derive a sharpened estimate of 
time-varying volatility that can be used in our multi-factor model. 



19 



Another related work is that of Ait-Sahaha et al. (2010), who describe cross-country 



returns using jump models based on the class of mutually exciting Hawkes processes. These 
generalizations of Poisson processes allow for a shock in one region of the world to increase 
the probability that future large shocks will be observed, both in the same region and in 
other regions. Market shocks are modeled as instantaneous jumps in the price path, rather 
than as conditionally Gaussian innovations with time- varying volatility. 

In the Hawkes-process model, the analogue of volatility is a latent vector of conditional 
Poisson jump intensities for each asset. These intensities are assumed to follow an vector- 
autoregressive-like process. Our model differs from this approach in two main ways: we 
model returns using stochastic volatility rather than a jump process; and we explicitly 
model the manner in which volatility shocks enter the cross-section of expected returns. 
By contrast, in a Hawkes-process model, there is no structural mechanism for changes in 
the conditional jump intensities (the analogue of volatility shocks) to affect the pattern 
of signs that one is likely to observe in the return residualsj^ There are also multivariate 
time-varying correlation GARCH in the mean models, such as the BEKK(1,1) model. We 
prefer to examine conditional volatility in a multivariate factor-model setting, much the 
same way as 



Ng et al. (1992). The main advantage then of our approach is being able to 
also account for directional volatility moves in returns, which we demonstrate in the next 
section to be crucial for quantifying contagion. 

4.2 Benchmarking of mean/ volatility models for the U.S. market 

We have described a model wherein heavy-tailed shocks to volatility may enter the cross- 
section of expected returns. We now present evidence that the severity of the shocks that 
arise during financial crises cannot be explained easily either in the Gaussian SV framework, 
or in the asymmetric GARCH-in-mean (GARCH- AM) framework. The explosive SV model, 
on the other hand, does quite well. 



^Regarding the issue of price-path discontinuities: when one marginalizes over the local scale factors At in 
our volatility process, one gets the increments of subordinated Brownian motion for at, which will also have 
discontinuous sample paths as the discretization becomes arbitrarily fine. Mathematically, the difference is 
between a finite-activity Poisson process superimposed on a diffusion process, and an infinite activity Levy 
process with a countable number of absolutely summable tiny jumps. For a formal discussion of this point, 



see Poison and Scott (20121 
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ESV versus GARCH-AM: EU Market, May 2010 
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Figure 2: Two estimates of daily volatility for the EU market during May 2010. 



Table 1: The ten largest residuals from GARCH-AM model fit to the U.S. market returns 
1963-2010 (top) and E.U. market returns 2000-2010 (bottom), compared with the cor- 
responding residuals from the ESV model. The model residuals are standardized by the 
estimated conditional volatility (thus 1 denotes a one-standard-deviation move). 



Date 


US Return {%) 


Model Residual 
GARCH-AM ESV 


July 7, 1986 


-2.86 


-4.92 


-2.37 


September 11, 1986 


-4.37 


-5.13 


-2.53 


October 19, 1987 


-17.16 


-7.76 


-4.91 


October 13, 1989 


-5.37 


-10.27 


-2.97 


November 15, 1991 


-3.43 


-6.84 


-2.36 


February 16, 1993 


-2.60 


-5.09 


-2.28 


February 4, 1994 


-2.29 


-5.46 


-2.20 


October 27, 1997 


-6.55 


-6.55 


-2.07 


January 4, 2000 


-3.99 


-6.05 


-2.37 


February 27, 2007 


-3.43 


-7.34 


-2.35 






Model Residual 


Date 


EU Return (%) 


GARCH-AM 


ESV 


January 4, 2001 


-4.08 


-3.56 


-2.32 


January 3, 2006 


3.52 


4.41 


2.38 


May 17, 2006 


-3.63 


-3.37 


-2.33 


February 27, 2007 


-5.60 


-6.72 


-2.70 


July 26, 2007 


-4.25 


-3.48 


-2.53 


September 4, 2008 


-5.23 


-4.19 


-2.76 


September 29, 2008 


-11.11 


-4.86 


-3.62 


April 27, 2010 


-5.18 


-4.12 


-2.69 


August 11, 2010 


-5.46 


-3.75 


-2.50 


October 19, 2010 


-5.37 


-4.43 


-2.61 
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Table 2: Results from the U.S. market model regressions. 



Model QQ X 100 oi X 100 Residuals 



Mean 


Volatility 


Est. [t-stat] 


Est. [t-stat] 


Skewness (p) 


Ex. Kurt, (p) 


Log-like 


1(a) 


GARCH-AM 

SV 

ESV 


3.0 [3.686] 
9.8 [9.630] 
8.7 [9.811] 


-0.2 [-0.190] 
-8.6 [-6.139] 
-6.8 [-5.693] 


-0.42 (0.000) 
-0.10 (0.003) 
-0.02 (0.565) 


2.62 (0.000) 
1.34 (0.000) 
-0.01 (0.806) 


-14029.97 
-12927.04 
-12833.40 


2(a) 


GARCH-AM 

SV 

ESV 


3.2 [1.802] 
16.4 [8.188] 
14.8 [8.365] 


-0.4 [-0.159] 
-16.5 [-6.303] 
-14.5 [-6.168] 


-0.42 (0.000) 
-0.11 (0.001) 
-0.06 (0.088) 


2.63 (0.000) 
1.13 (0.000) 
0.11 (0.015) 


-14029.97 
-12935.19 
-12853.19 


1(b) 


GARCH-AM 




2.4 [2.861] 


-0.42 (0.000) 


2.60 (0.000) 


-14037.04 


2(b) 


GARCH-AM 




5.6 [6.119] 


-0.41 (0.000) 


2.59 (0.000) 


-14033.45 


3 


GARCH-AM 


2.9 [4.703] 




-0.42 (0.000) 


2.62 (0.000) 


-14029.98 



The lag-1 asymmetric GARCH model (Glosten et al. , 1993) holds that today's volatil 



ity is a deterministic function of previous market shocks and volatility states. IVIoreover, 
volatiliy is assumed to respond asymmetrically to positive and negative shocks: 



Co + Cie?-i + C2eLi(l.*-i<o) + C^<yti , 



inducing returns to exhibit the so-called "leverage effect." When volatility appears in 
the model for expected returns, as here, we have an asymmetric GARCH-in-mean model 
(GARCH-AM). 

We fit these five mean models cross with the three volatility models (GARCH-AIVI, 
Gaussian SV, ESV) to the daily excess returns on the European and US market portfo- 
lios, considered independently. The maximum- likelihood estimate of the GARCH model 
was computed using gradient descent. To fit the SV models, we used a sequential IVIonte 
Carlo algorithm called particle learning (PL). For an explanation of the particle- learning 



algorithm, see Carvalho et al. (2010a), along with Lopes and Tsay (2011). Each SV model 
has three unknown parameters associated with the volatility process: the long-term mean 
of volatility (ctq), the persistence of volatility (0), and the scale of the volatility shocks 
(r^). In traditional empirical-Bayes fashion, we searched a grid of these values to find the 
combination that maximized the log-likelihood of the data for different model considered. 
For the ESV model, we fixed the degrees-of-freedom parameter to v = 2. 
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The data set used for estimating the U.S. market model runs from July 1, 1963 to 
October 29, 2010; comprises 11916 daily returns; and is freely available from Kenneth 
French's websitej^ The excess return on the US market is the value-weighted return on 
all NYSE, AMEX, and NASDAQ stocks, minus the daily return on a one-month Treasury 
bill. Likewise, the excess return on the Euro-wide market is the return on the Vanguard 
ETF corresponding to the pan-European index from Morgan Stanley Capital International, 
minus the daily return on a one-month Treasury bill. The data set used for estimation runs 
from March 11, 2005 to October 29, 2011; and is freely available from Yahoo! Finance, using 
the ticker symbol VGK. 

Table [T] and Figure [2] help to convey a sense of the differences that can arise between 
the ESV model and the asymmetric GARCH-in-mean model (G ARCH- AM), another state- 
of-the-art volatility estimator. (See Section [3] for details.) Both of these models produce 
an expected value fit and an estimated standard deviation at for the daily excess return, 
denoted xt- But these estimates for the state variables can be strikingly different. For 
example. Figure [2] plots the estimates of daily volatility (a^) from both models for May 
of 2010, when the Greek sovereign- debt crisis reached a climax and the assessment of the 
two models temporarily diverged. The main difference arose on the 10th of May, when 
EU finance ministers announced the creation of a 500-billion-euro loan package intended to 
stop Greece's debt troubles from spreading to other Euro-zone economies. The E.U. stock 
market promptly rose by 9.5%, the fourth-largest daily gain since the adoption of the euro. 
In response, the ESV estimate of volatility nearly tripled, and remained high for several 
days; over the same period, the GARCH-AM estimate remained roughly constant. 

Neither behavior, of course, is automatically more sensible for the particular events in 
question. But the residuals {xt — fJ't)/f^t provide a natural way to quantify the fit of the 
two models across an entire time series. By this standard, the ESV model behaves much 
more sensibly. Table [T] lists the 10 largest residuals for the GARCH-AM model for each of 
two data sets: the excess returns on the U.S. market portfolio from July 1963 to October 
2010, and the excess returns on the E.U. market portfolio from August 2000 to October 
2010. If returns are conditionally Gaussian given fit and at, there is a 50% chance of seeing 

^ http : //mba. tuck. dartmouth. edu/pages/f acuity /ken. french/Data_Libraxy/f-f .factors .html 
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a single move of five (or more) standard deviations in a period of 1000 years. If we believe 
the G ARCH- AM estimates for at, the U.S. market saw 9 such moves in 47 years. As the 
table also shows, these daily residuals are still large under the ESV model, but not absurdly 
so. 

This snapshot communicates the essence of the argument: when markets undergo drastic 
shocks, our estimate of volatility must change rapidly, especially if we hope to quantify the 
role of volatility risk in the cross-section of daily returns. 

Table [2] shows the results of a more formal model assessment, which reveals that the 
explosive stochastic volatility model, version 1(a), has the largest log-likelihood. (Models 
lb, 2b, and 3 were not competitive; we include their GARCH estimates for comparison.) 
We also use two other tests of model fit that assess the third and fourth moments of the 
normalized residuals, which should correspond to those of a standard N(0, 1) distribution. 
First, we applied D'Agostino's test for skewness, where the null hypothesis is that residuals 
have zero skewness. Second, we applied Anscombe's test, where the null hypothesis is that 
the residuals have no excess kurtosis, to detect tails that are systematically different from 
the standard normal distribution. Table [2] also reports the empirical skewness and kurtosis 
of the standardized residuals for each of the models, along with the two-sided p-values for 
D'Agostino's test and Anscombe's test. Both of these tests strongly contra-indicate the 
GARCH-AM and Gaussian SV models, but not the explosive stochastic-volatility model. 



5 Quantifying contagion in the EU sovereign debt crisis 

Using Model 1(a) in conjuction with an explosive stochastic volatility model, we estimated 
filtered means and volatility states for both the US and EU market risk factors. We used 



these results to construct the predictors in the four- factor model implied by Equation (10). 
The first two factors are returns on market-baskets, while the other two are our correspond 
to our filtered explosive- volatility factors: 



E{ya\EU, US) = /3f + /3f ^^f ^ + 1^'^^ + if^^. 



t ! 
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where yit is the return on index i, S^^ is the volatihty shock to the US market, and Sf'-^ 
is the excess volatihty shock to the European market. The excess shock is defined as the 
residual after regressing the EU volatility shock upon the US volatility shock. This measures 
the part of the EU volatility shock that cannot be explained by shared dependence; this 
is necessary to avoid marked collinearity, since the US volatility shock strongly predicts 
the EU volatility shock. In constructing our volatility risk factors, we use filtered volatility 
estimates, rather than smoothed estimates (which investors could not, in principle, know 
contemporaneously) . 

We first fit this model assuming static values of all factor loadings. These parameters 
are crucial in that they measure the market co-movements that can be explained by cross- 
section dependence upon volatility shocks. Recall, moreover, from Section 3 that there is a 
direct correspondence between the 7 parameters and the covariance matrix between return 
and volatility shocks. By estimating 7, we are reconstructing this covariance matrix. 

Our estimates of these parameters are summarized in Table |3j (In this table and others, 
"Euro" refers to the daily return on an exchange-traded funded that holds Euros in cash.) 
These results demonstrate excess correlation in the residuals after controlling for the first 
two factors (that is, relative to asset-pricing model that includes global and regional market 
risk). Notice the stark difference in magnitude of the 7^ coefficients when the ESV model 
is used to estimate factor volatilities, as opposed to the GARCH-AM model. Our findings 
suggest that there is a real, substantively meaningful effect of volatility shocks in the cross- 
section of returns — but that this effect is subtle, and can easily be masked unless one 
pays very careful attention to estimating volatility states themselves. A comparison of the 
magnitude of the coefficient estimates is particularly instructive here. Even though the 
GARCH-AM volatility model leads to many statistically significant loadings, the practical 
significance of these effects is largely muted, often by an order of magnitude or more. 

The main empirical utility of our model, however, is in studying the time-varying be- 
haviour of the volatility loadings Jj^'^jjf^, particularly during crisis periods. To explore 
this effect, we introduce use dummy variables to measure period- by-period variation in the 
loadings on the volatility shocks. We include dummies for the two most notable "EU- 
centric" crisis periods of May 2010 and August 2011. 
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EU Market Daily Returns and Estimated Volatility: 2000-08-01 to 2010-10-29 
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Table 3: The results from fitting the four- factor model incorporating contagion effects 
due to differential transmission of explosive volatility shocks. Loading estimates are listed 
with their t-statistics in brackets. Countries are listed in ascending order of their loadings 
on the EU residual volatility factor. We include estimates of model parameters under 
both the GARCH-AM and ESV models, to demonstrate how our conclusions regarding the 
importance of volatility in the cross-section only emerge under the superior volatility model 
(as measured by log- likelihood from the previous section). "Euro" refers to the daily return 
on an exchange-traded funded that holds Euros in cash. 



Explosive SV 






/ 












DEU 


0.125 [ 


4.30] 


0.903 


[38.16] 


0.444 [ 


5.44] 


-0.229 [ 


-2.79] 


Euro 


-0.443 [ 


-19.7] 


0.524 


[28.57] 


0.050 [ 


0.78] 


-0.17 [ 


-2.68] 


UK 


0.200 [ 


6.80] 


0.810 


[33.91] 


0.141 [ 


1.71] 


-0.180 [ 


-2.18] 


SWE 


0.177 [ 


3.88] 


1.084 


[29.22] 


-0.247 [ 


-1.93] 


-0.221 [ 


-1.71] 


BEL 


0.119 [ 


2.95] 


0.824 


[25.10] 


-0.290 [ 


-2.56] 


-0.162 [ 


-1.42] 


CHE 


0.057 [ 


1.93] 


0.719 


[29.78] 


-0.186 [ 


-2.23] 


0.011 [ 


0.13] 


ERA 


0.062 [ 


2.92] 


0.972 


[55.83] 


0.058 [ 


0.97] 


0.045 [ 


0.75] 


ESP 


0.003 [ 


0.08] 


1.047 


[37.52] 


0.018 [ 


0.19] 


0.097 [ 


1.00] 


NLD 


0.132 [ 


4.76] 


0.879 


[39.00] 


-0.236 [ 


-3.03] 


0.195 [ 


2.49] 


ITA 


0.023 [ 


0.70] 


0.994 


[37.48] 


-0.006 [ 


-0.07] 


0.25 [ 


2.72] 


GARCH-AM 






/ 








if 




DEU 


0.1 


[4.1] 


0.92 


[44.4] 





[-1.6] 





[-3] 


Euro 


-0.44 [ 


-21.8] 


0.52 


[32.6] 


-0.01 


[-3.1] 


-0.02 


[-2.5] 


UK 


0.23 


[8.8] 


0.77 


[37.2] 





[-1.1] 





[-0.9] 


SWE 


0.19 


[4.8] 


1.06 


[32.7] 


0.02 


[3.6] 


-0.01 


[-1.5] 


BEL 


0.09 


[2.5] 


0.84 


[30.4] 


-0.01 


[-1.3] 


-0.02 


[2.1] 


CHE 


0.06 


[2] 


0.71 


[33] 


0.01 


[2.2] 


-0.01 


[-1] 


ERA 


0.04 


[2] 


1 


[65.4] 





[1.5] 


-0.03 


[-5.9] 


ESP 


-0.07 


[-2.4] 


1.11 


[42.7] 


0.01 


[2.4] 


-0.01 


[-0.8] 


NLD 


0.13 


[4.9] 


0.88 


[44.3] 





[0] 


-0.01 


[-0.8] 


ITA 


-0.05 


[-1.5] 


1.07 


[44.2] 


-0.01 


[-2] 


-0.02 


[-4] 
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For comparative purposes, we will also include the loading estimates for the financial 
crisis of 2008. Late summer and fall of that year provided the height of the financial crisis. 
On Monday, September 15, 2008 the Dow Jones plunged 504 points with the US authorities 
trying to put a rescue package together for AIG for a $20 billion lifeline. On Tuesday 
September 30, Asian markets were the first to react to the shock that the $700 billion Wall 
Street bailout had failed to pass through Congress. On Wednesday, October 15, 2008, the 
FTSE suffered its fifth biggest fall in history, closing down 7.16% at 4079.5 (a 315 point 
fall); the Dow Jones dropped by 7.8%. The following day, the self-exciting nature of these 
shocks led (on Thursday, October 16) to Japan's Nikkei suffering its worst fall since 1987, 
and to the FTSEIOO slumping again by a further 218 points (to 3861). 

These events form a useful "control population," in that they are global rather than 
EU-centric, and thus we would expect to see different patterns of loadings on the volatility 
shocks during these two periods. To quantify the need for dynamic structure, we have 
also provided a p-value for the partial F-test corresponding to a test of the time-varying 
contagion effect model, versus no time-varying effect. 

The results from our basic pricing model are provided in Table [3j Clearly different 
countries load upon aggregate volatility shocks in different, economically interpretable ways 
during different crises; thus the factor loadings cannot be explained simply by nonlinear 
dependence of returns on volatility. For example, Germany has significant negative loadings 
on volatility shocks during the major crisis periods of May 2010 and August 2011. These two 
contagion terms have individually significant t-statistics, and the F-test strongly supports 
their inclusion in the model. 

An interesting comparison is with Spain and Italy. These two indices have significantly 
significant positive loadings on volatility shocks during 2010, and significantly negative 
loadings during August 2011. In the May 2010 period, Germany had to play a large role in 
bailing out Greece, raising the probability of a bailout for other debt-stricken nations. But 
in August 2011, the prospects for a satisfactory resolution to the Greek situation looked 
significantly worse, and market fears about Spain and Italy intensified. 

For France, there are a number of differences between the periods of May 2010 and Au- 
gust 2011. In particular, the major concerns over the French banks Societe Generale, BNP 
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Paribas, and Credit Agricole weren't present in 2010; all of these banks faced downgrades 
by ratings agencies in 2011, due to their massive exposure to Greek bonds. 

Also instructive is the pattern of loadings exhibited by Switzerland (CHF). During both 
the 2008 and 2011 crisis periods, the Swiss Franc was seen as a safe haven for investors — 
particularly in 2011, when fears gathered about the possible dissolution of the euro. (The 
franc went from 0.77 euros in early April to 0.97 euros on August 10, 2011; no such spike was 
observed during May of 2010.) The corresponding estimates during these periods suggest a 
noticeable "flight to quality" effect. Finally, Sweden — which has been relatively unaffected 
by the EU debt crisis — was the only country where the volatility shocks were consistently 
irrelevant. 

6 Conclusion 

The recent European debt crisis has brought to focus the inter-dependence of market move- 
ments. The increasingly explosive behavior of returns and volatility in periods of market 
turmoil has heightened the need for statistical models capable of describing contagion ef- 
fects, such as those eminating from of the Greek debt crisis. These effects are hard to 
measure, with the challenge arising from the facts that excess correlation can be caused by 
co-movemcnts in an aggregate risk factor, and that conditional factor means and volatilities 
cannot be observed directly. 

This leads us to consider dynamic models of the expected return/risk trade-off, a rela- 
tionship that is fundamental to any asset pricing model and to models which can uncover 
time-varying correlations. To uncover the effects of explosive volatility in explaining con- 
tagion, we propose a four factor model, fit to data from major European equity markets 
during the period 2008-2011. In this model, loadings on the volatility risk factors have nat- 
ural interpretations in terms of time-varying patterns of correlation in a mutually exciting 
stochastic volatility model (Section 3), whereby shocks to returns and shocks to volatility 
are correlated cross-sectionally. We construct the volatility risk factors themselves using a 
heavy-tailed explosive stochastic volatility model, which is comprehensively benchmarked 
against the popular GARCH-AM framework. 
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Our findings are two- fold. First, traditional volatility models are not sufficiently sensitive 
to extreme market movements, and are thus inefficient at estimating conditional market 
volatility. Secondly, after controlling for global and regional market integration, we still 
find evidence for contagion in the European debt crisis. We have gone beyond merely 
demonstrating the presence of excess correlation, however. Indeed, we have shown that 
some of this excess correlation can be explained by postulating volatility effects in the cross 



section of returns, paralleling the construction of Ang et al. (2006) for domestic equities. 
But these effects only emerge when the explosive SV model is used; they are largely muted 
when the demonstrably inferior G ARCH- AM model is used instead. This foregrounds the 
need for efficient, state of the art volatility estimation techniques in studying presumptive 
periods of financial contagion. 

One deficit of our model is that it is purely statistical nature, and is therefore silent on 
the underlying macroeconomic causes of such events. However, we feel that this data period 
and our multi-factor approach provides a benchmark for future empirical finance studies 
of contagion, particularly given the readily interpretable nature of our findings (e.g. Table 
|4]). Moreover, as the literature review in Section 2 explains, there are many theoretical 
reasons why the study of aggregate market volatility, and the transmission of volatility 
shocks, should play an important role in understanding contagion. These theories — which 
encompass a myriad of potential causes, from behavioral to macroeconomic to financial — are 
buttressed by the fact that the explanatory power of volatility shocks in the cross section 
appears to be strongest during periods of crisis. 
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